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AERODYNAMICS OF A RECTANGULAR WING OP INFINITE ASPECT 
RATIO AT HIGH ANGLES OF ATTACK 
AND SUPERSONIC SPEEDS 

By John C. Martin and Frank S. Malvestuto, Jr. 

SUMMARY 


Pertijrbatlon of the flow over a two-dimensional flat plate at finite 
angles of attack is used to obtain a first-order evaluation of the damping 
in roll, the lift and moment due to an increment in angle of attack, and 
the lift and moment due to a steady pitching velocity for a rectangular 
wing of infinite aspect ratio at supersonic speeds. Approximate expres- 
sions are derived for the lift and moment due to a constant vertical 
acceleration . 

The results are valid for the ranges of Mach n\miber and angle of 
attack for which the flow behind the shock is supersonic. The analysis 
is based on the eq.tiations for rotational flow, so that the change of 
entropy is taken into accoiant. 

Design charts are presented which permit rapid estimations to be 
made of the aerodynamic derivatives for a given Mach number and a given 
angle of attack. 

The res'ults for the Infinite-aspect-ratio wing are used to make 
estimates of a number of the aerodynamic derivatives for rectangular wings 
at finite angles of attack. 


INTRODUCTION 


The development of the linearized theory of supersonic flow has 
permitted a first-order evaluation of a nmber of aerodynamic stability 
derivatives for a variety of plan forms at an angle of attack of 0°. 
Second-order theories similar to the one introduced by Biisemann (ref. l) 
and extended by Van Dyke (refs. 2 and 5) have been used to obtain second- 
order evaluations of a number of the aerodynamic derivatives for a few 
simple airfoils (refs, k to 6) at supersonic speeds. The second -order 
theories predict no variation in the aerodynamic derivatives with angle 
of attack. There is a need, however, for values of aerodynamic derivatives 
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at angles of attack beyond the validity of the second-order theories 
developed in references 4 to 6. At the present time, little published 
information can be found on the aerodynamic derivatives at finite angles 
of attack and supersonic speeds. 

The only publisihed papers associated with these aerodynamic deriva- 
tives that have come to the authors ’ attention are the analyses of Ivey 
(ref. 7), Car-rier (refs. 8 and 9), and Chu (ref. 10). The analyses by 
Carrier and Chu make use of the linear perturbation theory for rotational 
flow (refs. 11 to 13 ) which allows first-order estimates to be made of 
the flow variables behind a strong shock attached to the leading edge of 
a two-dimensional, wedge or within the region bounded by the lower surface 
of an airfoil at finite angles of attack and the strong shock from the 
leading edge of the airfoil. 


The present paper contains a first-order evaliaation of a number of 
aerodynamic derivatives for a rectangular wing of Infinite aspect ratio 
at finite angles of attack, based upon the linear perturbation theory for 
rotational flow. This analysis. Including the development of the line- 
arized rotational-flow equations, was performed independently of previous 
analyses, an attempt being made to present a completely unified treatment 
leading directly to the evaluation of aerodynamic stability derivatives . 
Wherever possible, similarity of results from the present and previous 
analyses are noted. 


The resTilts are valid for the ranges of Mach number and angle of 
attack for which the flow behind the shock from the leading edge is 
supersonic. A first-order evaluation of the following aerodynamic deriva- 
tives is made: the lift-curve slope Cj^, the rate of change of pitching 


moment with angle of attack Cpi^, the damping in roll 


Cjpt/ the lift due 


to constant pitching the moment produced by a constant rate of 

pitch Cioq^. 


Simple approximate relations for 


CniQ_, Cjpi, Cj^, and 

are derived. These approximate relations yield results which are in good 
agreement with the exact first-order values, except at angles of attack 
near the angle where the flow behind the shock is sonic. In addition, 
approximate expressions are determined for Ci^ and Cm^, the lift and 


pitching moment due to a constant vertical acceleration. It should be 
noted that, although the shock-expansion theozy can be used to calculate 
and (see ref. 7) with relatively little effort, the use of 


this theory to evaluate the remaining derivatives becomes difficult, if 
not impossible. The methods used herein yield the first-order evaluation 
of the aerodynamic derivatives for the airfoil considered with relatively 
little effort. 
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A series of design charts presented herein permits rapid estimations 
of the aerodynamic derivatives to be made for a given Mach number and a 
given angle of attack. 

The resiilts for the wing of infin ite aspect ratio are used to make 
estimates of a number of the aerodynamic derivatives for rectangular wings 
at finite- angles of attack. 


SYMBOLS 


A 


aspect ratio 


K- 


a = 1 + 


■n 


Mi^r(r - 1) 


B= \ 

- 

1 

Bo = ' 


- 1 

Bi = 


- 1 

B2 = 


- 1 


b 

bi,b2,b5^b4 

c 

°1 


°2 

°v 

G(x,y,z) = 0 


wing span 
constants 

chord 

velocity of soiznd behind shock 

velocity of so-und in flow over upper surface of airfoil 
specific heat at constant volume 
equation of pertxirbed shock surface 


g 




acceleration due to gravity 

arbitrary functions of (x - B 2 _z) and (x + Bj_z ) , respectively 
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J 

K 

% 

Kii 

M 

Mo 

Mq_ 

M2 

m 

P 

pP' 

pq 

pa 

pa 

P 

P' 


Po 


Pi 


enthalpy in free stream 

enthalpy in flow behind shock from leading edge 

unit vectors in x-, y-, and z-direction, respectively 

mechanical equivalent of heat 

constant 

parameter defined by equation (47) 
parameter defined by equation (52) 
parameter defined by equation (92) 

Mach nxmiber 

free-stream Mach number 

Mach number behind shock 

Mach mmiber of flow over upper surface 

slope of shock (see fig. 5) 
pressure coefficient 

pressure coefficient due to a constant rate of roll 
pressure coefficient due to a constant rate of pitch 
pressure coefficient due to an Increment in angle of attack 

pressure coefficient due to a constant vertical acceleration 

pressure 

rate of roll 

free-stream press-ure 

pressure behind shock in unperturbed flow 
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P 2 pressiare in unperturbed flow over upper surface of airfoil 

6p-j^ first-order perturbation in pressure of flow behind shock 


q. rate of pitch 

S entropy 

8Si first-order perturbation in entropy 


s distance along shock from leading edge of airfoil in 

xz-plane 

T temperature 

t time 


u,v,w perturbation velocity components in x-, y-, and z-directlon, 

respectively 


q = 

Vo 

Von 



Vi 

V2 

Vi' 


iu + jv + kw 

free-stream velocity 

component of free-stream velocity that is nonnal to the 
shock 

component of free-stream velocity that is tangential to 
the shock 

velocity of unperturbed flow behind shock 

velocity of unperturbed flow over upper siarface of airfoil 

velocity vector of perturbed flow behind shook 



normal component of velocity in pertxarbed flow behind shock 



tangential component of velocity in pertiirbed flow behind 
shock 

ccHEponent of in xy-plane 

component of in xy-plane 
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Wma-y maximum value of W 

W = (Vi + u)l + vj + wk 

vertical velocity of airfoil associated with Ai 

Xcg distance from leMing edge of airfoil to center-of-gravity 

location 

x,y,z rectangular coordinates 

ct, angle of attack 

Oq fixed vBlije of a 


a = 


5g 


motion of the wing corresponding to a constant vertical 
acceleration} sometimes referred to as a plunging motion 


r 

9 

Sy 

P 

Po 

Pi 

P2 

8pi 

$ 

0 


ratio of specific heat at constant pressure to specific 
heat at constant volume (l.4O0 for all calculations) 

angle between free -stream direction and unperturbed shock 
profile (see fig. 5) 


density 

free-stream density 

density in unperturbed flow behind shock 

density in unperturbed flow over upper surface of airfoil 
first-order perturbation in density behind shock 

scalar potential fimctlon defined by equation (97) 
scalar potential function 

scalar potential fxmctlon associated with rolling 
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scalar potential function associated with a steady pitching 
velocity 

scalar potential, function associated with perturbation 
in a 

scalar potential function associated with a constant 
vertical acceleration 


V= i— +j— — 
OX oy oz 


^ = iBi^ L. + j L. + k ^ 
Sx ^y 3z 


"N 


nonnal-force coefficient. 


Force 


1 2 

— P^Vq X Plan-form' area 




pitching-moment coefficient. 


Pitching moment 
1 2 

— c X Plan-form area 


a 


roUliJg -moment coefficient about stability axis, 

Rolling moment 

1 ? 

oPo^o ° ^ Plan-form area 


Ct. = 


^ / 


cos Oq 


o?=ao 






\ 


cos 


vOC 


% 


o£r=ao 
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cos Oq 


Note that ^L(X^ *^Lq. defined herein are the component deriv- 

atives noirmad to the free stream. 




Sa /cB=a^ 



Cm^ - 



^l-g 



Wherever the variables x, y, z, and t are vused as STjbscrlpts, 
a differentiation process with respect to these variables is indicated. 


PRELIMINARY REMARKS 


The airfoil considered in this paper is a rectangular wing of 
infinite aspect ratio at finite angles of attack (fig. 1) . The airfoil 
is taken to be thin so that the thickness effects can be neglected in 
the first-order evaluation of the aerodynamic derivatives formulated 
herein. The aerodynamic derivatives are obtained by finding the first- 
order perturbation in the flow over a two-dimensional flat plate at a 
finite angle of attack (fig. 1(a)). The reader is assumed to be familiar 
with the shock-expansion theory of two-dimensional supersonic flows. The 
stability axes used in the analysis are shown in figure l(b) . 
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OLq, *-1%, <ierlvatives are determined by considering the 

effect of an infinitesimal increment in a. The Cjpi derivative is 

determined by analyzing the effect of an Infinitesimal constant rate of 
roll about the x stability axis. 


The CLq^ and 0^^ derivatives are determined by considering a con- 
stant infinitesimal rate of pitch about the y stability axis. Since the 
type of motion analyzed in finding and is often misunderstood. 


this motion is discussed in detail. A constant rate of pitch is associated 
with a constant rate of rotation about the axis of pitch while the angle 
of attack with respect to the free stream remains constant. Viewed from 
a point fixed with reference to the undisturbed air, the wing is flying 
in a circle of radius Vp/q. with a constant angular velocity q and with 
a constant angle of attack (fig. 2). Note that the flow associated with 
this motion is steady. Since the rates of pitch considered are very small, 
the raditis of the circle is very large and the basic flow for the airfoil 
considered herein can be taken to be the steady flow over a two-dimensional 
flat plate. 


The approximate expressions for the Ol^ and derivatives are 

determined by ass-umlng a constant infinitesimal acceleration in the 
z-dlrection of the stability axes. 


In general, the boundary conditions on the airfoil surface are com- 
plicated by the presence of the shock and expansion fan from the leaxilng 
edge of the wing. As an example, take the case of a small Increment in 
the angle of attack, which can be considered as the result of a small 
constant vertical velocity Wg, stq)erimposed on the original flow. In 

the stability-axes system (which is fixed relative to the airfoil), the 
free-stream direction has changed by the angle tsi. (fig. 3) • Since the 
shock and the expansion fan at the leading edge change the direction of 
the flow so that it becomes parallel to the airfoil surface, the normal 
component of the pertxirbation velocity on the surface must be directly 
proportional to Ax times the unperturbed flow parallel to the surfane. 
The normal component of the pert\irbation velocities is related to the 
pngTp. of attack only through the unperturbed flow parallel to the surface. 
The boundary conditions on the surface due to other types of motions can 
be determined in a similar manner. 


PERTURBED FLOW OVER UPPER SURFACE 


The basic unperturbed flow over the upper surface is miform and 
Irrotational. The flow begins with the termination of the expansion fan 
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around the leading edge of the airfoil. InaBmuch as the flow is irrota- 
tional and the small disttirbances produced oh the upper surface will not 
Interact with shock waves upstream of the trailing edge of the wing, these 
distinrbances can be calculated by the tiBe of a potential function. 

The coordinate axes used in the analysis of the flow over the upper 
surface are indicated in figure k. The potential function must satisfy 
the following partial differential eq^uations; 

For steady flows, 

+ J^yy + 5^zz = ° 

For unsteady flows, 

-B2^0xx + ?5yy + ^ i^xt - ?^tt = 0 

The subscript 2 refers to the conditions on the upper surface of the 
airfoil. 

The boundary conditions associated with the various motions are: 

Upstream of the end of the expansion fan from the leading edge of the 
airfoil, 


^ = 0 


For an Increment in a, 



-Ax V2 


For a constant rate of pitch, 

2^0 = - ^Cg) 

For a constant rate of roll. 
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For a constant vertical acceleration, 

&=0 -o>V2t 

The potential functions associated with the various motions are; 
For an increment in a. 




Ax. V2(x - B22 ) 

B2 


For a constant rate of pitch, 





2 

- B 2 z) _ q.Xcg(x - Egz) 
2B2 


where the axis of pitch is located a distance X(,g fran the leading edge 
of the wing. 

For a constant rate of roll, 

^p« ^ pV(x - Bgz) 

B2 

For a constant vertical acceleration. 


A 

B2 


M2^^ M-2^ o t \ 

— + z + V 2 t(x - B 2 z) 

2B2^^ ^ 


The pressure coefficient "based on the free -stream conditions is 
given by 
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or 


P = 


2Mg/P2p2 fe , 

Mo V^o^Vo VgVo/ 


(2b) 


It follows that the perturbation presstire on the u^pper surface of the 
airfoil resulting from the varioaos motions is: 

For an increment in a, 

P” . . ^^2^ 2 (3a) 

BaMo^Po 

For a constant rate of pitch. 




qc_ 

2Vo 


4M2 

jV^2 X - Xcg 

B2Mo\ 

1 PoPo 


For a constant rate of roll, 


pp ' _ slz ^ 

“ ' 2Vo PoPo 

For a constant vertical acceleration, 

^ l ^2P2^ x ^2^^'*^ 


(5b) 


(5c) 


(5d) 
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STEADY TWO-DIMEWSIOKAL PERTURBED FLOW 
BEHIND M INCLINED SHOCK 


The perturbed flow over the lower surface of the airfoil due to 
steady pitching and a small Increment in a is a special case of the 
first-order, steady, two-dimensional pertiirbation of the flow behind a 
two-dimensional inclined shock. This flow may not be iirotational; 
therefore, the analysis is based on the equations for rotational flow. 

As stated in the introduction, e3q>ressions for the evaluation of first- 
order pe 2 ?t^lrbatlons of flow variables behind a strong shock have been 
derived in a number of papers. The develoiment and application herein 
of the linearized equations of rotational flow are an effort to present 
a complete unified treatment appropriate for the evaluation of aerodynamic 
derivatives. The coordinate axes and some other data used in the analysis 
are Indicated in figure 5- 

The equations for rotational flow are (p. 202 of ref. l^l-) ; 

The Eioler equations of motion. 



(V X W) X ¥ = - ivp 
P 


(M 


The equation of continuity. 


V.pW = 0 (5) 

The entropy eqiiation. 


- ( 6 ) 

The equation for conservation of enthalpy. 


7 P _ ¥§ax (Y) 

7 - 1 p 2 2 

The preceding equations are presented in their three-dimensional forms 
because these are needed in the analysis of the rolling motion. Equa- 
tions (4) to (y) represent six equations for five variables (p, p, and 
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the three components of W) since the entropy S is determined from the 
boundary conditions. It follows that there are only five Independent 
relations contained in these equations. The six equations are retained, 
however, because the sixth equation is tised to shorten the analysis. 

The flow behind the shock can be expressed as; 


W = i(Vi + u) + Jv + kw (8a) 

P = Pi + 8pi (8b) 

p = P;j^ + bp]^ (8c) 

S = (8d) 


where u, v, and w are perturbations in the velocity and SPi# 

and 5Si are the first-order pertiirbations in density, in pressure, 

and in entropy, respectively. The subscript , 1 refers to conditions 
on the lower surface of the airfoil. 

Substituting equations (8) into equations (it) to (7) and retaining 
only the first-order terms yields (these operations are given in the 
appendix) : 


Vu + 



^ 

8y/ \8x 8z / 


1 

ViPi 


VBpi 


(9) 


8x 8y 8z Pi 8x 
5Pi 6p^ BS^ 

"pI “ ~pI" 


( 10 ) 


( 11 ) 
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spi spx y ■ 1 ViPi sa?i 
Pi ' Pi ■ 7 Pi "" Ti 

Inasmuch as the entropy Is constant along streamlines, the change of 
entropy 6Si is a function of z only. The initial values of 6Si 

at the shock are part of the boundary conditions on the shock. 

When the perturbation flow is two-dimensional, v and all quantities 
operated on by — are zero, and eqmtlons (9) and (lO) reduce to 


Vn + 




V8pj_ 


(15) 


^ dz Pi dx 


(14) 


The partial differential equations for u and w can be obtained 
as follows: Since 6Si is not a function of x, the partial derivatives 

of equations (ll) and (12) with respect to x are 


and 


i = 0 (15) 

Pi Bx Pi Bx 


1 ^5Pl 1 ^5pi ^ y _ i ViPi 

Pi Bx Pi ox 7 p^ Bx 


(16) 
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By eliminating and from equations (l4) to (l6) and 

by lasing the relation 



the follQ-wing equation can be obtained; 


^ 1 ^=^ (17) 

ox 02 

Equation (I7) is one relation between u and w. A second relation can 
be obtained as follows; The pearbiarbation pressure c£in be eliminated 
from the z-component of eqaaations (ll) and (15) to yield 


^1 Pi ^1 


(18) 


The pearbiarbation pressure can also be eliminated from equations (ll) 
and (12) to give 


?£l = _ ^1^1 ^ _ ^1 

Pi 7Pi (7 _ l)c-v 

Equations (18) and (19) can be combined to yield 

^ ^ ^ P^s ^^1 

8x 8z V2 ^P]^Cy( 7 - 1) 8z 


(19) 


( 20 ) 


Equation (20) is the second relation between 


u and w. 
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The following partial differential equations for u and w can 
he obtained from equations (I 7 ) and ( 20 ): 


B^2 5^u _ _ 


Pi 




3x' 


2 ^2^ VqPlC-yCy - 1) 3 z 2 


( 21 ) 


.2 


w 


w 


Bi^ — - 


= 0 


(22) 


The general solutions of equations (21) and (22) can be written as 


u 


g^^^x-B^z) gg^x+B^z) p^_ 8 S^(z) 

Bi ViPl(7 - l)cv 


Bn 


( 23 ) 


w = gj^(x~B±zj 


g 2 (x+Biz) 


(2k) 


From equations (U) and (l2), the perturbation pressure is related 
to u by 


Pn 6 St 

5Pl = - ViPiU (25) 

(7 - l)cv 

Also from equations (U) and (12), the perturbation density is related 
to u by 


_ %^Pl^ Pi ^1 

■ Tfl '(/-Dct 

Equations (23) to ( 26 ) are the solutions of equations (U) to (l4). The 
functions gj_, gg^ and bSq must be determined from the boimdary con- 
ditions on the surface of the body and on the shock. 
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The boundary condition on the surface is given hy 


(w)z=o = 'Vif(x) (27) 

where f (x) is a known function of x determined hy the type of dis- 
turbance being investigated. 

It should be noted that although the shock surface is free to move 
and deform, the boundary conditions are not satisfied on the actual sur- 
face but along the shock profile of the unperturbed flow. This approxi- 
mation is v a l i d within the bounds of the linearization process used 
herein. 

The boundary conditions on the shock must be determined from the 
basic shock eq^uations. These equations are (from pp. 97 98 of 

ref. l4) : 


Po^on - Pi^l^i' 

( 28 ) 

Po'^on^ + Po = Pl('^ln‘)^ + Pi 

( 29 ) 

PoVot = 

(50) 

From equations (28) and ( 30 ) it follows that 

(51) 

V = ^lt‘ 

( 32 ) 


Since (from p. 99 of ref. l4) 


y Pq 
7 - 1 gPo 


Jho = 
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and 


Jli]_ = 


y Pi 

7 - 1 gPl 


equation ( 51 ) can be expressed as 


_2L_ + 1 Vo 2 

7 - 1 Po 2 


7 - 1 Pi 2 V / 


(55) 


The normal and tangential velocity conrponents can be written as 
(see fig. 6) ; 


Von = ■'^o sin 0 = Vi’sln(0 - a) 

Vo^ = Vo cos 0 Vi^' = Vi'cos(0 - a) 

Substituting these four relations Into eqiiatlons (28), (29), (52), 
anl (55) yields 


PqVq sin 0 = p^V^'sln(0 - a) 


p^Vo^sln^© + Po = P3^^V3^’^^sln2(e - a) + 

Vo cos 0 = V2^'cos(0 - a) 

7 ^ Vo^sln^0 _ y Pi (vi‘) sln2(9 - g) 

7-1 Po 2 7-1 Pi 2 


(5i^) 

(55) 

(56) 

(57) 


The first-order variations across the shock wave can be determined 
by taking the total derivative of the preceding equations while considering 
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Po, Po^ Vq as constants. The differentials dpj_, d0, da, 

and dV^' can he considered as first-order variations of the variables. 

Thus, from the four preceding equations, the differentials are related 
by the following equations; 

PoVo cos 0 d0 - V]_'sln(0 - a) dp^_ - Pj^ sln(0 - a) dV^' - 
PlVi'cos(0 - a)(d0 - da) = 0 (38) 


2 

2pQVo^sln 0 cos 0 d0 - sin^(0 - a) dpj - 2 p^V 3 _’sin ^(0 - a) dVj_' - 

2p^^V^'^ sin(0 - a)cos(0 - a) (d0 - da) - dp^ = 0 (39) 

-Vo sin 0 d0 - cos(0 - a) dVi’ + Vi*sin(0 - a) (d0 - da) = 0 (4o) 

Vorsin 0 cos 0 d0 ^ — + —21— ^ ^ - V, 'sin2(0 - a) dV. ' - 

7 - 1 Pi 7 - 1 Pi Pi 

^Vi'^^sin(0 - a)cos(0 - a) (d0 - da) = 0 (4l) 

5i:om figure 6 it can be seen that 

dVi' = u 



and 
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By ^ing use of these relations and 
to (4l) can be expressed as 


equations (54) to ( 57 ), 


equations ( 58 ) 


(Po - Pi) Vo cos 0 d9 - ^ u sin 0 + w cos 0 
\ / Vl Vi 


sin 0 dpT = 0 

Pi 1 


(42) 




o 'o 


2V. 


^ sin20 _£_ 0 gQg 0 _j_ 


2 

Vq^ Sin 20 dpj_ + dp^ = 0 


(^3) 


° Vo Sin 0 d0 + ^ u cos 0 + ^ ^ w sin 9 = 0 (44) 


Pi 


Pl^l 



sin 0 cos 0 d0 - 


^ u sln2e + PoVof 

^1/ PqVi 


w sin 0 cos 0 - 


y ^Pi ^ r Pi<^i 

7 - 1 0l 7-1 " ° (Its) 


and, "«=-antlala 

i:ised. This relation is ' relation between u and w will be 
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where 



Note that for zero angle of attack 



Bo 


Values of Kj are gi-ven in table I for various angles of attack and 
Mach numbers. 

Equation (^7) is one of the boundary conditions on the shock. It 
remains to determine the values of 6Si in the flow behind the shook. 

Since is constant along streamlines (and hence is a function of z 

only) the values of 8Si are determined by the conditions on the down- 
stream side of the shock surface. 

The ejqjression for Si/cy is (from p. 202 of ref. ik) 


Si 

— = logg Pi - 7 loge Pi 


m 


Taking the total differential of equation (it-8) yields 


dS 1 5S 1 dpi 

Cy Cy Pi 


6pi 

7 

Pi 


(^9) 


This equation is true for points located directly behind the shocks thus 
for points on the downstream side of the shock, equation (49) becomes 
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5Si dp2^ dp2_ 

Cy 


(50) 


Expressing dpi and dp^^ In terms of u by the iise of equa- 
tions (42) to (45) gives for equation ( 50 ): 




u 

Vl, 


z=mx 


(51) 


where 




Mo 


2 Pq/P1 

Pl\po 


2 

) (r 



^ ^ Bin^e 
Pi Pi 




(52) 

Values of Kj are presented in table I for various angles of attack and 

Mach n-umbers. Equation ( 5 I) conclx^ies the expressions for the boundary 
conditions. The other expressions for the boundary conditions are equa- 
tions ( 27 ) and (46) . 

The boundary condition on the lower surface of an airfoil due to a 
small increment in the angle of attack is (in the coordinates shown in 
fig- 7) 


(w)2=0 = Ax Vi = gi(x) + g£(x) 


( 55 ) 


It follows from equation (53) that 


w = da Vi 
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and from equation (46) that 


u = Kj Ai 

The perturbation pressure is (frcm eqs- (26) and (54)) 


(5iJ-) 


5Pi = -PiV3_ 


Kj Ac Vt + 


Pi 5Si(z) 


(55) 


Since (from eqs. ( 5 I) and (54)) 


6 S 1 

= KiiKi A. 


eqioation (55) can be written as 


5Pi = -PiVi^Ki At 


K. 


1 + 


n 


Mr 7(7 - 1) 


(56) 


This expression for the pressure has previously been derived by Chu 

, Kii 


(ref. 10). The quantity Kj 1 + 


1 


Mi^7(7 - 1 ) 


I 


equals -tan A in equa- 


tion (43) of Chu's paper. The increment in the pressure coefficient 
based on the free-stream conditions is 


P_ 

Ax 





Mi^r(7 - 1) 


(57) 


Figtnce 8 presents the variation of 


p/d 


dP/Ax 


with angle of attack for various 
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Pi Pi 

Mach numbers. The values of Mt , — , — . and 0 were taken from ref- 

^ Po Po 

erence 15 and from unpublished calculations. 

The effect of the change in entropy in the perturbed flow can be 
evaluated by setting Kjj equal to zero in equation (57)* The change 

in entropy is retained in the boundary condition on the shock which 
relates the velocity components u and w. Figure 9 presents a com- 
parison of values of ^/Ax with and without the change in entropy 
Included in the perturbed flow. 

An approximation to the pressure due to an increment in a can be 
obtained by assuming that the presence of the shock in the iierturbed flow 
can be neglected and by assuming that the perturbed-flow velocity com- 
ponents can be expressed by derivatives of a potential function. This 
potential function is similar to the function associated with the flow 
over the upper surface and is based on the velocity Under these 

asstanptlons, the velocity component u can be expressed as 


u 


Ax Vx 

"bI~ 


( 58 ) 


The pressure coefficient associated with equation ( 58 ) and based on the 
free -stream conditions is 


p ^ 2MiSl 

^ BlMo^o 


( 59 ) 


Figure 10 presents a con^iarlson between exact and approximate values of 
the increment in the pressure coefficient on the lower surface of an air- 
foil for an increment in a at various Mach mmabers. In view of the 
assumptions involved in the approximation, figure 10 indicates that equa- 
tion ( 59 ) is a very good approximation to the pressure coefficient except 
in a very small angle— of— attack range near the an^le of attack for which 

Ml = 1. 


PERTURBED FLOW OYER LOWER SURFACE DUE TO STEADY PITCHING 


The boundary condition on the lower surface of a pitching airfoil 
is (in the coordinate system shown in fig. 7) 
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(w)z=0 = q.(x - Xcg) = gi(x) + g2(x) (60) 

If g^(x-Biz) g2(x+Biz) for z 0 are as Binned to "be of the form 

gl(x-Biz) = gbiCx - B^z) + b2 

g2(x+Biz) = qb^Cx + Biz) + bi^. 

where b^^, b2, bj, and b^ are constants, then u and w are given 

by (from eqs. (23), (24), and (51)) 

w = q(bi + b5)x - q(bi - b^)Biz + b2 + bij. (6l) 

(62) 

zpsmx 


u = - 


q(bj_ - h^) , t )2 - ■*=4 

— X + q(bi + bj) z 


K 


n 






Mi^rCr - 1) 


u 


From equation (6o)/ 

q^x - Xcg) = q(bi + b?)x + b 2 + b^ 
Equating powers of x yields 

b 2 + b^ = -qxcg 


bi + bj = 1 

On the shock (from eqs. (46), (6l), and. (62)), 

q(bi - bjjmz q(b2 + ^3)2 b2 - bj,. ^ 

B;^a a B 2 _a 


q(bi + b 5 )mz - q(bi - bj^Biz + b 2 + b^ 


Ki 


(63) 

(64) 


(65) 
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where 


a = 1 + 


K- 


•n 


- 1) 

Equating powers of z In equation ( 65 ) yields 


1 + KjBq_a 1 - KjBQ_a 

^2 + — = 0 


Bqa 


Bia 


(66) 


1 + KtBt a 1 — KxB-i a 

Eia (“ - ®l)”l - Bia (° + = ° 


{ 67 ) 


Equations ( 63 ), (64), (66), and ( 67 ) axe four equations containing 
the four unknown constants bi, b 2 > bj, and bi(.. Solving for the 

unknown constants yields 


2^m - KjBiSaj 


bo = — 




J2 2~~\ ' ^1^1®’ 


‘) 


— Bji^ ^ ^1 + KjB^^a^ 
^ 2^m - KjBi^aj 


A 

2 V 


bi^ = - ^{1 + KjBia) 


The velocity components w and u are 
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w = 


m - 


- Rjina^B^^ 


- qxcg 


( 68 ) 


u = - 


m 


- KjBi^a L: 


^1 - X - - KjB3_2a^ z - qx^gK^-i 


ra + 


qKijKi 


m2 - Bt2 


Mi^yCr - 1) 


^cg 


( 69 ) 


m 




The pressure on the surface is (from eqs. (25) and (69)) 


m - Kja 

Spi = pjViq — X + XcgKja 

1 - KiBiTaa 


(TO) 


This expression for the pressiire can eilso be readily extracted from 
Chu's analysis (ref. lO) by substituting in equation (72) of that paper 
qx^Vi for f'(x), the condition requiring the flow to be tangential to 

the airfoil surface. 


The pressvire coefficient on the lower surface of the airfoil is 


p = 3£_ 5 + Kra Si (71) 

2Vo Mo yPoPo 1 - K^Bi^nia ° 

Figure 11 presents the chordwise pressxire di^ribution on the lower sur- 
face of an airfoil for vario\as angles of attack at Mq = 2.00 where 

the axis of pitch is located at the midchord point. 

The effect of neglecting the change in entropy in the perturbed flow 
can be obtained from equation (7I) by replacing a by unity. The change 
in entropy is retained in the boundary condition on the shock which relates 
the velocity components u and w. Thus, the pressure coefficient on 
the lower surface of a pitching airfoil with the change in entropy neg- 
lected is 
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P = S2_ ^1 /PlPl / ^ ~ X jj. 

2^0 ^ f PoPo\“ “ ° 

Figure 12 presents the chordwlse variations of the 
including and neglecting the change in entropy for 

angles of attack with the axis of pitch located at 
This f Igvire indicates that the change in entropy has a strong effect on 
the pressure for the higher angles of attack. 

An approximation to the pressure due to pitching can he obtained 
by making the same assiHirptlons as were, made for the approximation to 
the pressure due to an increment in a. This procedure will lead to the 
following relation for u: 


( 72 ) 


pressure coefficient 
Mq = 2.00, for two 

the mldchord point. 


- ^cg) 


(73) 


The pressure coefficient associated with equation (73) based on the 
free-stream conditions is 


P w 


qc ^1 I PlPl /x 

^^ypoPo\'= 



(7^^) 


Note that equation (7^) can also be obtained by replacing Kj by - l/B^ 

in equation (72). Figure 15 presents exact (from eq. (71)) approxi- 
mate (from eq. (7^)) values of the chordwlse pressxire on the lower surface 
of a pitc hing airfoil for two angles of attack at Mq = 2.00 with the 
axis of pitch located at the midchord point. This figure Indicates that 
equation (jif) is a good approximation to the pressure distribution 
for oq = 9 . 7 ° and 20.7°. 


PERTUECBED FLOW OVER LOWER SURFACE DUE TO ROLLING 


The perturbed flow over the lower surface of an airfoil due to rolling 
can be determined from equations ( 9 ) to (l2). The coordinate axes used 
are shown in figure 7 where the y-axis is directed into the plane of the 
page. 
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Tbe boijndary condition on the surface is 


(w)z=o = p‘y 


( 75 ) 


The boundary conditions on the shock are determined from the eq^ua- 
tions expressing the conservation of mass, momentum, and energy. The 
continuity of mass flow req^uires that (eq. ( 28 )) 


Po^On Pl^ln 


(t6) 


It can be seen in figure l4 that 


Vq^ = Vq sin 9 cos(dX) 


= Vq' Sin(0 - a) cos(dA) 

To the first order in dA the preceding equations become 


Vo^ = Vo sin e (77a) 

Vln' = Vq' sin(9 - a) (77b) 

Thus, to the first order in dA equation ( 76 )^ becomes 

PqVq sin 0 = PqVq' Bin(0 - a) ( 78 ) 


This equation is the same as the corresponding relation for the two- 
dimensional case. 

The conservation of momentum in the direction normal to the shock 
requires that (eq. ( 29 )) 
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i^oni eq,uations (TT), the preceding equation can be written (to the first 
order In dX) as 

PoVo^sln^e + Po = p3^(Vi')^Bln2(0 - a) + p^^ (80) 

This Is the same aa the corresponding relation for the two-dimensional 
case. 

The conse 2 rvation of the tangential momentum requires that (eq. (30)) 

PoVoa^ot = Pl%’''lt' 

It follows from equation ( 76 ) that 

Tot = Tit' 

Figure l^(a) shows that the tangential components in the xz-plane are 

(Vot)x2 “ ® 

Thus (eq. ( 36 ) ), 

Vo cos 0 = Vi*cob (0 - a) (81j-) 

This is the same as the corresponding relation for the two-dimensional 
case. 

From the preceding resrilts, the energy law can be expressed as 
(eq. (37)) 


2 

y Po Yo^sin^O _ y Pi ^ (Vi* ) sin^(0 - a) 

7-1 Po 2 7-1 p^^ 2 


(83) 


4l 


This is the same as the corresponding relation for two-dimensional flows. 
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The equation for the conservation of mass flow, the equation for 
the conservation of momentum in the direction^ normal to the shock, the 
energy equation, and the equation of the tangential velocity components 
in the xz -plane are the same as the corresponding equations for the two- 
dimensional case. It follows from the analysis of these relations given 
for the two-dimensional case that two of the boundary conditions on the 
shock are given by equations (^4^) and ( 51 )* 

Figure l4(b) shows that the tangential velocity components in the 
y-direction are (to the first order in dX) 

(Vot) = Vq sin e dX (86a) 

y 


= Vi sin(e - a) dX - V (86b) 

Thus , from equation ( 82) , 

Vq sin 0 dX = Vj sin(0 - a) dX - v 
or 


V = “1^0 9 - sln(9 - a)JdX ( 87 ) 

The slope of the shock surface in the y-dlrectlon dX is given by 


dX=^ (88) 

where G(x,y,z) = 0 is the equation of the perturbed shock surface. 
This surface G can be expressed as 
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where s is the distance along the ehock surface in the xz -plane measured 
from the leading edge of the airfoil. From eqixations ( 87 ) , ( 88 ) , and ( 89 ) , 
the velocity component v can be expressed as 


V = 



sin 0 


- sln(9 



( 90 ) 


From equations {h2) to (^5) > "the velocity component u and d0 are 
related by 


d0 = Kttt 


V- 


^/z=mx 


(91) 


where 


K 


1 . Mo2 50 Po 

I Pi Pi ! 


Ill 


iW^ 

\Po 


cot% 


1 H ) cot 0 

^ Po; 


1 + J: _ Mq 2 Pol sin20 (l - 2 ^ 

Po Pi V Pi; 


(92) 


Values of Kjjj are given in table I for various angles of attack and 
Mach numbers. Substituting eqxiation ( 9 I) into equation ( 90 ) yields 


cLs (95) 

z=mx 

Equations (46), (51)^ (95) are the boundary conditions on the shock 

surface . 

The following relation can be obtained in a manner similar to that 
tised in deriving eq-uatlon (21) : 


V = - 


Vq sin 9 - Vq sln(0 - a) 


n^iii 

IQ 

/ 

^ ^1 

Jq \^yj 


Sx By Bz 


( 94 ) 


4 
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The components of equation (9) are 


3u _ 

1 


Sx 

VlPl 

Sx 

bv _ 

1 

S&Pl 

Sx 

Vi 

Sy 


1 

S5p^ 

Sx 

ViPi 



Eliminating 8p-j_ from equations (95) yields 

^2 v2 

O U O V 



dx2 

_ 


Sx Sz 

dx^ 

a^v _ 



^ Sx Sy 


( 95 a) 


(95b) 


(95c) 


These equations can be expressed by the vector equation 
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Since the vector is irrotational, it can be described in terms 

of a scalar potential function where 



(97) 


From equation (9 ^) ) it follows that 


Vj,. 0 (98) 

ox 

Substituting eq\aation (97) into equation ( 98 ) yields the following 
partial differential equation for $: 


+ «yy + = 0 (99) 

Thus, the scalar potential function (j must satisfy the three-dimensional 
wave equation. 

The boundary condition on the surface for $ is (from eq. (75)) 

The boundary conditions on the shock surface are expressed in terms of u, 
V, and w and are not easily expressed as conditions on 

A solirtion of equation (99) which satisfies the condition given by 
equation (lOO) is 


«D = Ky 

where K is a constant. If the conditions on the shock can be satisfied, 
the partial derivatives of the velocity components with respect to x 
are given by 


^ _ du 
3x ^x 



36 


NACA TK 3^21 


by bx 


bil> _ bv _ Q 
bz bx 

Thus, the velocity components u and w are independent of x 
he written as 

u(x,y,z) = u(mz,y,z) 

■w(x,y,z) = w(mz,y,z) 


The ho\mdary condition on the surface req^uires that 


(v)z=o = p'y 


Thus, the velocity component v must be of the form 


w = p'y + n(z) 


where Jl(z) is an unknown function of z which is zero when z 
The velocity con^nent u is given (from eqs. (46), (lOla), and 
tiy 


u = Kijp’y + n(z)j 

Since 


bv _ 
bx 


then by integration 


V — (v)jj— JJJ 2 + K(x — mz) 


and can 

(101a) 

(101b) 

( 102 ) 

is zero. 
( 102 )) 

(103) 
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From equations (93) ^ (lOlb), and (103)^ 


(■v)x=mz = [vo sin 8 - Bin(e - a<)J KjjjKj^mS + 1 p'l 


Thus, 


= |vq sin 0 - V 2 _ sin(9 - a) + 1 p’z + K(x - mz) 


( 104 ) 


Substituting eqxiations (102), (103)^ and (104) into equation (94) yields 
the .relation 


Thus, 


5v _ dfl(z) _ Q 
3z dz 


n(z) = 0 


because 0(0) = 0. 

The velocity con^jonents u and w are 

u = Kjp’y 

w = p'y 

The change in entropy is given (from eqs. ( 51 ) and ( 105 )) by 


( 105 ) 

(106) 


SSi _ KjjKjp'y 

<=v 


( 107 ) 


The pressure is (from eqs. (25)^ (105)? and ( 107 )) 
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8pi = -ViPj^Kjp'y 


pi%i%p'y 

(7 - l)Vi 


(108) 


From equations ( 95 t>)# aJ^d (108) and the relation 


K 

8 x 


the constant K is found to he given hy 


K = Kjap’ 


Thus, equation (l 04 ) becomes 


V = |vq sin 0 - Vj_ sin(0 - a)^ ^ \jxS + 1 p’z + K^;ap'(x - mz) 

^ (109) 

The pressure coefficient based on the free-stream conditions is given by 


pi^ 

p'y Mo ]j PoPo 

2Vo 

Figure 15 presents the variation of with angle of attack for 

p’y/2Yo 

various Mach numbers. 

The effect of the change in entropy in the perturbed flow can be 
evaluated by setting ■ Kj j eqiial to zero in equation (UO). The change 

in entropy is retained in the boimdary condition on the shock which 
relates the velocity components u and w. An examination of equa- 
tions (57) and (no) indicates that the percentage changes in the pres- 
sure coefficients which result from neglecting the change in entropy are 
the same for a small Increment in a and for a small rate of roll. 


K- 


1 + 


II 


Mi^/Cy - 1) 


(no) 
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p 

Thus, the effect on — - — ■, of neglecting the change in entropy can he 

P y/2Vo 

ohtained from, the values of /!p/Ax given in figure 9- 

An approximation to the pressure due to rolling can be obtained by 
neglecting the presence of the shock and assuming that the perturbed- 
flow velocity coioponents can be expressed by derivatives of a potential 
function. This potential function is based on the velocity a.n(^ is 

similar to the potential function associated with the flow over the i 5 >per 
svirface. Under these assimptions, the velocity component u can be 
expressed as 


u«-£lz 

Bl 


( 111 ) 


The pressinre coefficient associated with equation (ill) and based on 
the free-stream conditions is 


(u^) 

pV BiMoyPoPo 

2Vo 

Figure l6 presents a comparison between exact and approximate values of 
the pressure coefficient on the lower surface of a rolling airfoil for 
various Mach numbers. This figure indicates that equation (112) is a 
good approximation to the pressure coefficient, except in a very small 
angle-of-attack range near the angle of attack where = 1. 


AEEODINAMIC DERIVATIVES FOR RECTANGULAR WING OF 
INFINZDE ASPECT RATIO 

The expressions for the pressure coefficients on the upper and lower 
surfaces of an airfoil d'ue to the various motions considered permit the 
calciilatlon of the aerodynamic derivatives associated with these motions. 

The derivative C]^ for an airfoil at an arbitrary angle of attack 
will be defined as 
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■'a 




cos Oq 


/ac=ao 

It follows from eqiiations (3a) and (57) that Cj^ is given by 


_2_< 

2 Pi 

1 

H 

1“ 

\ 

2 

. !!2_£2L 

Mo^ 

^■1 — 
Po 

w 

Ml2y(7 - 1) . 

Po 

0 


(113) 


Note that the center of pressure remains at the midchord point. Figure 17 
presents the variation of GLo, with angle of attack for various Mach 

numbers . 

The derivatives and for an airfoil are (from eqs. (5b) 

and (71)) 


n 2 

Mo 


M , iPl / m - Kj-a ^ ^I^cg ^ PgPg /^l ^ 
^'^PoPoll - KiaBiSj c B2 ' c / 


COB Oq 

(llij.) 


Cm^ 


4M- 


1 PlPl 


Mo VPoPo 


m 


~ ^ 1 _ ^cg\ ^ Kjaxcg /i _ ^cg^ 


1 - KiaBi2m\3 2c 


c \2 c 


^^2 /P2P2 


3MoB2V^oPo 


1 _ 5 ^£8 + 5 
c 




(115) 


Figures I8 and 19 present the variations of Cl^ and respectively, 

with angle of attack for various Mach numbers when the axis of pitch is 
located at the midchord point. 
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The derivative C^, of an airfoil is (from eqs. (3c) and (UO)) 






?lPl 


'PoPo 


1 + 


K- 


■II 


Mi^r(7 - 1) 


- ^ ^ (n6) 

BaVPoPo 


Figure 20 presents the variation of C 2 1 'vd.th angle of attack for 
varioi;is Manh numbers. 


APPROXIMATIONS TO AERODYNAMIC DERIVATIVES FOR WING OF 
INFINITE ASPECT RATIO 


The approximate expressions for the pressure coefficients on the 
lo'wer surface of an airfoil permit the deri'vation of simple approximate 
expressions for the aerodynamic coefficients associated with the motions 
considered previously. 

An approximate expression for Cj^ is (from eq.s. (3a) and (59)) 

CCr 




2 Al^ ^ ^ 5^2^ 5^ 

Mo^lfl Po ^ ^0/ 


cos Oq 


(117) 


Figure 21 presents a comparison between exact and approximate values 
of Ct for variolas Mach numbers. This figure Indicates that equa- 

tion ( 117 ) yields a good approximation to the exact linearized value 
of given by equation (H3)> except in a very small angle -of -at tack 

range near the angle of attack for which M]_ = 1. 


Approximate expressions for Cj,^ and are (from eqs. (3b) 

and (74)) 
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4* 





°o 


(118) 




4 1% ^iPl 
SMoWyPoPo 


+ 




1 



(119) 


Note that the approximation for C^q (eq. (ll8)) yieldE a value of zero 
when the axis of pitch is located at the midchord point. 


Figure 22 presents a comparison between exact and approximate Values 
of Cjj^ over a range of angles of attack for an axis of pitch located 

at the midchord point. This figure indicates that equation (ll9) is a 
good approximation to eqioation ( 115 )^ except near the angle of attack 
for which Mi = 1, when the axis of pitch is. located at the midchord 
point . 


An approximate expression for 



is (from eqs. (5c) and (ll2)) 


CZp' 


5Mo\Bi 


/ PlPl ^ Mg , 
PoPo 


Eg’ 


P2^' 

PoPo/ 


cos 


( 120 ) 


Figure 25 presents a con5)arlson between exact and approximate values 

of C? t for angles of attack to the point where Mi = 1, for various 
"P • 

Mach numbers. This figure indicates that equation (120) is a good approxi- 
mation to equation (ll6), except near the angle of attack where Mj_ = 1. 


The results of the approximate expressions for the pressures and 
aerodynamic coefficients indicate that these expressions are in good 
agreement with the exact values except for angles of attack in a small 
range near the angle of attack where = 1. An approximate expression 

for the pressure on the lower surfane of an airfoil .jrfiich has a constant 
vertical acceleration can be obtained by neglecting the presence of the 
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shock aoad by •using a potential fimction based on the velocity 
This potential f'unction is given by 



Mi^x^ 

+ 

2Bi2 2 


+ V^t^x - 



( 121 ) 


The approximate pressure coefficient based on the free-stream conditions 
is, for t = 0, 


P ^ ^ % / PlPl X 

MopoPo ° 
2Vo 


( 122 ) 


Approximate expressions for Cj^ 
and (122)) 


and are (from eqs. (3d) 


2 A^l felPl ^ fepP2\ 

Cj^. sa _■£. -±- L , A - - i: cos Oq 

^\bi^|/PoPo 


^ MBi5pcf>o 

Figure 2Jf presents the variation of the approximate values of 


or — — -with angle of attack for varioiis Mach numbers. 

1 ^cg 

3 ' 2c 


( 125 ) 


( 124 ) 

cos Oq 


The 


sum of the aerodynamic coefficients and Cm^ partly 


determines the damping of longitudinal oscillations of aircraft and is 
directly proportional to the aerodynamic damping of slowly oscillating 
airfoils (see appendix B of ref. 6); for these reasons the sum Cm + Cm^ 
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is considered separately. Figures 25 to 28 present the variation 
of (calculated from, the approximate relations, eqs. (ll9) 

and (124)) with angle of attack for various Mach numbers and center-of- 
gravity locations. These figures indicate that, for an axis of rotation 
located between zero and the midchord point, Qmq + is positive 

(destabilizing) for some part, or all, of the angle -of -attack range for 
which the results apply. Figijre 28 indicates that, for an axis of rota- 
tion located at the three-quarter -chord point, Cnig^ + is negative 

(stabilizing) for all angles of attack for which the results apply. 


ESTIMATES . OP AEROD'XNAMIC DERIVATIVES FOR RECTANGULAR 
WINGS OP FINITE ASPECT RATIO 


The development of exact and approximate expressions for the aero- 
dynamic derivatives of the wing with infinite aspect ratio permits 
estimations of the aerodynamic derivatives for rectangular wings. These 
estimations can be based on either the exact or the approximate expres- 
sions for the inf inite-aspect -ratio wing. The approximate expressions 
will be vised because they are much simpler and they yield results which 
axe within a few percent of the exact values most of the angle -of - 
attack range for which the expressions will be used. 

The pressure within the tip regions can be approximated by assuming 
that the chordwlse and spanwise variation of the pressure in this region 
is the same as the pressure variation determined by linear theory for 
small angles of surface deviation. The Mach number associated with the 
values from the 31near theory is taken to be the Mach number of the flow 
over the surface in the region which is unaffected by the tips. Thus, 
the Mach number associated with the linear-theory values for the flow 
over the upper surface is M2> and the Mach number associated with the 

flow over the lower surface is M]_. With these assumptions, the following 

approximate expressions for the aerodynamic derivatives based on the 
approximate expressions for the wing of infinite aspect ratio and the 
exact linear -theory expressions (from ref. l6) can be written as 


Ei/i. _ 1 \ ^ £ 2 / 0 ^ _ 1 \ 

,^ 0 / Po\pi 2ABi^/ \f^/ Po\®2 2AB2V 




cos Uq 


(125) 
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n ' « J fe 1 ^11 2 

^ 1 W Po 3ABi 


i + 2!Eefi__L_ 

c \ 2 AB 2 


_ 1 + 2 _ 


P2 1 2 


2ABiy Po B2 5AB2 


ci,^ ^ i . 4 ^ssA . + 

MoVPqPo®! 5ABi c I 2ABii 


M 2 /P 2 P 2 1 P^2 " 2 1, Xcg/, 1 


_£ /1£L£. £ - ■ _ k :iHS. 1 ±_ 

MoV^oPq ®2 3AB2 ° \ 


cos 



1 

5 - 8ab^ 


5ABi 


--#8. 





1 

3 - 8AB2 

2B2 

3AB2 


+ :i£B 8 - JL 

C \ AB 2 


^28) (8 - JL 


C « _ fe^/2 _ ^ _1 ^ 1 \ ^ 

^p' 2 BiMo\/PoPo\5 ABi 5a2bi2 i2a5bi5/ 


^2 rMkfi. i_ + 1 ^ 1 

AB2 5A^B2^ 12A-^^ 
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/£lPl _1_(_2 + ^ '*' 


^\l^oPoB2p\ 3ABi 



2P2 1 I 2 + ^ '*' 


5AB2 


cos Oo 


(150) 


Ml 1 L 2^ 2 + Bi^ ^ ^ + 2Bi^ j . 

^ Mol/poPo Bi3ly3 " c ■ 2ABi c 3^^! j 

^ -L .4 . 2 ^ ^ •*•■^2^ + ^cg (151) 

Mo^PoPoBg^l^ c 2AB2 ° 5AB2 j 


The expressions for the aerodynamic derivatives based on the lin- 
earized theory which were xised in eq,mtlons (125) "to (l5l) limited 
to cases where the parameter AB is eqxial to or greater than 1. Thus, 
equations (l25) to (131) are limited to angles of attack and Mach nimi- 
bers where ABi is equal to or greater than !• For angles of attack 

and Mach numbers that cause ABi to lie between 1 and Z, the tips affect 

over one-half the wing area on the lower surface; thus the values cal- 
culated from equations (125) to (13I) may differ considerably from the 
true val\;ies. Figures 29 to 39 present the variations of the estimated 
aerodynamic coefficients of rectangular wings with angle of attack for 
aspect ratios of 2 and i|-. In these figiares the dashed portions of the 
curves indicate the region where the parameter ABq lies between 1 and 2. 


DISCUSSION OF EESULTS 


If the change in entropy in the perturbed flow is neglected and the 
change in entropy in the boundary conditions on the pei*turbation velocity 
components on the shock surface is retained, a large error results for 
Mach numbers of 3 snd 4. The results also indicate that the aerodynamic 
derivatives for the wing of infinite aspect ratio can be approximated 
very well by simple expressions over most of the angle-of -attack range 
for which the theory is valid. 
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The results for the aerodynamic derivatives seem to indicate that 
there is a rough analogy between increasing the angle of attack at a 
given Mach number and decreasing the Mach number at zero angle of attack. 
For example, the approximate expression for Cl^ given by equation (II 7 ) 

can be written as 






Po 



B=Bi 

ce=0 


+ 



Po 2 b=B2 

Of=0 


As the angle of attack increases, the first term on the right-hand side 
of the preceding equation increases, and the second term on the right- 
hand side decreases. Inasmuch as q ^•PProS'Ches infinity as the 

Mach ntimber approaches 1, this type of variation will dominate the varia- 
tion of equation (II 7 ) as the angle of attack approaches the angle for 
which Ml = 1. An examination of the other approximate expressions for 

the aerodynamic derivatives shows similar analogies between increasing 
the angle of attack at a given Mach number and decreasing the Mach num- 
ber at zero angle of attack. 

The calculated veilues of the aerodynamic derivatives show very rapid 
changes with angle of attack near the angle of attack where M]_ = 1. It 

is expected that, provided the lower surface of the airfoil is not flat, 
the thickness of the airfoil may tend to modify the rapid changes with 
angle of attack becaiise the thickness of the airfoil woiild change the 
local flow over the lower STirface. 

The variations of the estimated aerodynamic coefficients of rectang\ilar 
wings with angle of attack presented in figures 29 to 59 indicate that the 
effect of aspect ratio on the variation of these coefficients with angle 
of attack is quite strong. These variations can be e3q>lained on the basis 
of the results from the wing of infinite aspect ratio at finite angles of 
attack and the results from rectangular wings at zero angle of attack 
(linearized theory) . For example, consider the variations of the esti- 
mated Cl^ with angle of attack for Mq = 4.0 as shown in figirre 29 • 

At the higher angles of attack, the changes in the pressiires, the den- 
sities, anrt the local Mach numbers in the basic flow cause the contribu- 
tion to the estimated of the lower surface to outweigh greatly the 

contribution to the estimated of the upper surface. The region 

Influenced by the tips on the lower surface Increases as the angle of 
attack increases j thus, the tips tend to decrease the estimated at 
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the higher angles of attack. This effect causes the curves of the esti- 
mated Cj^ at the higher angles of attack and the higher Mach numbers 

to have slopes considerably lower than those of the curves for the wing 
of infinite aspect ratio. 


CONCLUDING RMABKB 


Perturbation of the flow over a two-dimensional flat plate at finite 
angles of attack is used to obtain a first-order evaluation of the damping 
in roll, the lift and pitching moment due to an increment in the angle of 
attack, and. the lift and pitching moment due to a steady pitching velocity 
for a rectangular wing of infinite aspect ratio at supersonic speeds. 
Approximations based on the results of the linearized theory and on the 
flow over a two -dimensional flat plate at finite angles of attack are 
derived and are shown to yield results which are in good agreement with 
the exact first-order theory. Approximate expressions are also derived 
for the lift and pitching moment due to a constant vertical acceleration. 
Estimations of the aerodynamic derivatives for rectangular wings of finite 
aspect ratio are derived by the combined use of the approximate relations 
for the Infinite-aspect-ratio wing at finite angles of attack and the 
restilts of the linearized theory for vanishingly small angles of attack. 

The same type of analysis -used herein can be applied to the swept- 
back wing of infinite aspect ratio with supersonic leading edges. The 
results for the aerodynamic derivatives of this wing could also be used 
to make estimations of the aerodynamic derivatives for a number of flnite- 
aspect-ratlo wings in certain Mach number ranges in a manner similar to 
the method tised in the present paper for the rectangular wing. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 25, 1955* 
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APPENDIX 

DETERMINAITION OF FIRST-ORDER-PERTURBAIEION EQUATIONS 
FOR FLOW BEHIND A TWD-DIMEIKIONAL SHOCK 


(eg. 


Consider first the vector form of the Euler equations of motion 

( 4 )): 


2 

V^ + (vxw) xW=.-^\;;^) (Al) 

The first teim of the preceding equation can be written, by the use of 
equation (8a), as 

V s£ , 

2 \ 2 22/ 


. i/vi ^ + u ^ + V ^ + w s»\ + 

^ 8x Sx 8x y 

j(vi ^ + u ^ + V ^ + w ^"\ + 

\ 8y By By By y 


k 



Bu 

Bz 


+ u 




V 



When only the first-order terms are retained, this expression becomes 



Vh 


(A2) 
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The second term of eqiiation (Al) can be expressed, by the use of 
equation (8a), as 


(V X W) X V = i 


w 


du 

dz 



V ^ + V §a) + 


(Vi + u)^ - fVi + u')^ - w ^ + w ^ 
''^x ' 'by by Sz 


k 


V 


bv 



When only the first-order terms are retained, this expression becomes 


(V X W) X W = i(0) + + kVife - 

v5x byj Szy 

The last tenn of equation (Al) cein be expressed as 


(A5) 




piV Pi 


• V6pj^ 


When only the first-order terms are retained, this expression becomes 


_ i Vt = - — V6pn (a4) 

P Pi 
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To the first order, eq.\iation (Al) can now be expressed, from equa- 
tions (A2) , (A5 ) , and (a 4) , as 


Vli + 



\Sx 



V6pi 

VlPi 


(A5) 


This is equation (9) in the body of the paper. 

Consider now the equation of continuity (eq. (5))* This equation 
can be expressed as 


V-pW = V*^Pi + ^^1 + ^ 


= 0 


(a6) 


When only the first-order terms are retained, the preceding equation can 
be expressed as (eq. (lO)) 


3u ^ Bv ^ Sw _ ^8pi 
hx Sy Pi ox 


(A7) 


The entropy equation can be expressed as 


p^ + 5p^ = ^pj^ + ^Pl)^ ® 


3j+SSi 
7 p ‘V 


or 


Pi + 5Pi = p/ + r 8Pi p/-^ ^ + . . . 


(A8) 


Pi - pi^ 


since (from eq. (6)) 
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Thus, to the first order in pressure, density, and entropy, equation (a 8) 
becomes equation (ll); 


Spj^ &p2_ 

P]_ ^ Px 


(A9) 


Now equation (A1) can be written, by the use of equations (8), as 


+ u2 ^ ^ 

■ ■ — + + — + — = 

7 - 1 2 2 2 2 

or as 

r / „ M . \ ."^1 . v2 , w2 W§ax 

^(Pl-sp,)-(i.— .. . .j. ^ 

(AlO) 


When only the terms \q) to and including the first-order are retained, 
equation (AlO) becomes 


r - 1 Pi r - 1 Pi 


7 ^1 

r - 1 Pn 


Pi 5pi ^ V- 


+ V^u = 




max 


(All) 


Since 

7 ^ ^1 _ ^'^ax 

7 - 1 Pi 2 2 


equation (All) can be reduced to equation (l2); 


^Pl _ ^Pl _ 7-1 u = Sk 
Pi " Pi ” Pi ®1 
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Figure Illustration of change in free-stream. direction viewed from 
the stability axes for a given increment in angle of attack. 
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Figure 6.- Velocity components associated with perturbed flow behind a 

two-dimensional shock. 



Figure 7*“ Coordinates used in analyzing flow over lower surface of an 

airfoil. 
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Annie nf attack rten 

Figure 8.- Variation of AP/6a with angle of attack on lower sirrface of 
an airfoil for various Mach numberB. 
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Figure 9.- Gangparlson of values of AP/Ai, on lover surface of an airfoil 
with and without change of entropy included in the perturbed flow. 
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Figure 10.“ CcanparlBon between exact and approximate valusB of AP/A 3 . on 
lower surface of an airfoil for various Mach numbers. 
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Figure 12.- Chordwise variation of pressure coefficient on lower surface 
of an- airfoil pitching about its midchord point, with and without change 
of entropy Included. Mq = 2.00j cxq = 9*7° 20.7°» 




Percent chord 

Figure 15.- Catrrparlson between exact and approximate values of chordwise 
pressure on lower surface of an airfoil pitching about its midchord 
point. Mq = 2.00; ttg = 9.7° and 20.7°. 






Figvtre 15 .- Vaxiatlon of pressure coefficient cm lower surface of a rolling 
airfoil with angle of attack for various Mach numbers. 
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Figure l6.~ Con 5 »arison between exact and approximate values of pressure 
coefficient on lower surface of a rolling airfoil for various Mach 
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(a) Mo = 1.25 to 2.50. 

Figure 17.- Variation of Cj^ with angle of attack for -varlcme Mach nmriberBc 
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Figure I9.- VariatlOD of ' 

"bera -when axis of pitch 
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Figure 25." CcsBaparison betweea exact aud approximate values of 

various Mach, nunibers. 



for 


VO 


MCA TN ^k21 




O 4 8 !2 16 20 24 2 8 32 36 40 


Angle of attack, deg 


(a) Mq = 1.50 to 4.00. 
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Figure 24.- Variation of approximate values of or := — 

cos Oq ► 1 _ ^cg 

5 ’ 2c 

angle of attack for various Mach numbers. 
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(b) Mq = 1.25 to 2.25. 
Plguxe 24.- Continued. 
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(a) « 1.25 to 1.75. 

Figure 25*- Variation of Cj^ + CW (calculated from approxljaiate eqs. (II9) 

and (l2i<-)) with angle of attack vhen center of gravity is located at 
leading edge . 
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Angle of attack, deg 

(b) Mq = 1.75 to 2.50. 
Figure 25.- Continued. 
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Figure 25.- Concluded. 
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Angle of attack, deg 


(a) Mq = 1.25 to 1.75. 

Figure 26.- Variation of 0^^ + Cm^ (calculated frcm approximate eq.a, (U9) 

and (124)) with angle of attack ^en center of gravity is located at 
quarter-cliord point . i 
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(b) Mq = 1.75 to 2 . 50 . 
Figure 26 .- Continued. 
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Angle of attack^ deg 

(b) Mq = 1.75 to 2.5(5. 


Figure 27.- Continued. 
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AngiG of attack, deg 

(a) Mq =! 1.25 to 1 . 75 » 

Figure 28.- Variation of (calculated frcm approximate eqa. (II9) 

and ( 124 )) with angle of attack vhen center of gravity is located at 
three-quarter-chord point. 






(b) Mq « 1.75 to 4.00. 

rigure 28.- Concluded, 









(a) Aspect ratio, 2.0, 

Figure 50 ■- Variation of estimated Cjj^ of rectangular wings with n ngl e 
of attack for various Mach numbers. Aspect ratio, 2.0 and 4.0. 
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Angle of attack^ deg 

(a) Aspect ratio, 2.0. 


Figure 51*- Variation of estimated 



of attack for various Mach numbers. 


of rectangular wings with angle 
Aspect ratio, 2.0 £ind k.O. 





(t) Aspect ratio, Iv.O. 
Figure 51.- Concluded. 
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Angle of attack, deg 

(a) Aspect ratio, 2.0, 

Figure 32.- Variation of estimated Cjjj^ of rectangular wings with angle 

of attack for various Mach numbers. Aspect ratio, 2.0 and 4.0j center 
of gravity located at quarter-chord point. 
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Figure 55 .- Variation of estimated of rectangular wings with angle 

of attack for various Mach, numbers. Aspect ratio, 2.0 to “j center 
of gravity located at mldchord point. 
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(a) Aspect ratio, 2.0. 


Figure 54.- Variation of estimated of rectangular wings with, angle 

of attack for various Mach numbers. Aspect ratio, 2.0 and 4.0. 
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("b) Aspect ratio, 4.0. 
Figure 5^.- Concluded. 
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(a) Aspect ratio, 2.0. 

Figure 55.- Variation of estimated of rectangular wings with angle 

of attack for various Mach numbers. Aspect ratio, 2.0 and if.O. 
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(b) Aspect ratio, 4.0. 
Figure 55 •- Concluded, 
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(a) Aspect ratio, 2.0. 


J'lgure 37 ... Variation of estimated Cm. 

of attack for various Mach numbers?^ 
o gravity located at midchord point. 
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(a) Aspect ratio, 2.0. 

Figure 58 .- Variation of estimated of rectangular wings for 

various Mach numbers. Aspect ratio, 2.0 and 4.0; center of gravity 
located at quarter-chord point. 
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(a) Aspect ratio, 2.0. 


Figure 59* - Variation of estimated Cjn^ + of rectajigular wings for 

various Mach numbers. Aspect ratio, 2.0 and 4.0; center of gravity 
located at midchord point. 





